Dynamics of a Bose-Einstein 
condensate in optical trap 

F.Kh. Abdullaev f , B.B. Baizakovf , and V.V. KonotopJ 
f Physical- Technical Institute, Uzbek Academy of Sciences, 
2-b, G. Mavlyanov str., 700084 Tashkent, Uzbekistan 
^Departamento de Fisica and Centro de Fisica da Materia Condensada, 
Universidade de Lisboa, Complexo Interdisciplinary 
Avenida Profesor Gama Pinto, 2, P-1649-003, Lisbon, Portugal 

February 1, 2008 



Abstract 

The dynamics of a 2D Bose-Einstein condensate in optical trap is stud- 
ied taking into consideration fluctuations of the far-off-resonance laser field 
intensity. The problem is described in the frame of the mean field Gross- 
Pitaevskii equation with randomly varying trap potential. An analytic 
approach based on the moments method has been employed to describe 
the noise induced evolution of the condensate properties. Stochastic para- 
metric resonance in oscillations of the condensate width is proved to exist. 
For the condensate with negative scattering length of atoms, it is shown 
that the noise can delay or even arrest the collapse. Analytical predictions 
are confirmed by numerical simulations of the underlying PDE and ODE 
models. 



1 Introduction 

Recently a Bose-Einstein condensate has been realized in all-optical far-off- 
resonance laser trap jjj . This important achievement opens up new perspectives 
in exploring the phenomenon. On the one hand, with the optical trap it is now 
possible to create more dense condensates, which is valuable for investigation 
of three body decay processes, to obtain different geometrical configurations in- 
cluding quasi-lD and 2D structures. On the other hand, in an optical trap atoms 
in arbitrary hyperfine states may be confined, therefore magnetic properties of 
atoms can be studied. 

However, alongside with the above significant advantages, optical traps have 
also some drawbacks, the most pertinent of which is related to fluctuations of 
the laser field intensity. The last circumstance introduces stochasticity into 
dynamical behavior of the condensate, which has to be taken into account in 
real situations. 
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In the mean field theory of BEC based on the Gross-Pitaevskii (GP) equa- 
tion, fluctuations of the laser field intensity can be regarded as modulations 
of the harmonic trap potential. The dynamics of BEC under regular varia- 
tions of trap parameters has been considered by many authors (see e.g. pi and 
references therein). From the experimental viewpoint the trap parameters are 
modulated in order to reveal the response spectrum of the condensate. Kagan 
et al., H have developed the scaling transformation of the order parameter, 
which in 2D leads to easily solvable ODE's giving results very close to exact 
numerical solution of the GP equation. The existence of parametric resonances 
in the condensate response to time dependent perturbations, on the ground of 
variational methods and numerical simulations, was reported in Ref. B. 

In the present paper we study the dynamics of a BEC under randomly 
varying optical trap potential, caused by fluctuations of the laser field intensity. 
The basic idea consists in employing of differential relations between several 
integral quantities of the GP equation. This approach, known as moments 
method, will be applied to analysis of the evolution of the atomic cloud's width 
driven by fluctuations of the laser field intensity. The issues of particular interest 
include the possibility of the stochastic parametric resonance in oscillations of 
the condensate width, and whether the noise induced expansion of the atomic 
cloud can delay or arrest the collapse of the BEC with attractive interaction 
between atoms. 



2 Basic equations and the Method of Analysis 

Starting point to further development will be the fact, that the dynamics of 
a nonlinear evolution equations (i.e. systems with infinite number of degrees 
of freedom) can be exactly reduced to the dynamics of a system with finite 
degrees of freedom ||. Namely, by introducing collective coordinates. This fact 
has been explored in Ref. j| in order to illustrate the existence of extended 
parametric resonances in a two-dimensional GP equation. An advantage of the 
model proposed in |J is that up to certain extent it allows exact solutions. 

We will concentrate on radially symmetric solutions (ip = ip(r, t)) to 2D GP 
equation, which describes the BEC having a cigar shape 

^ = ~\^/'t + (1 + £(<))r ^ + X ^ + ^ (1) 

The dimensionless variables of this equation are related to actual quan- 
tities through the following transformations: r — > r/a,ho, t — > tu>/2, ip — > 
(87rA f |a|/a/ lo ) 1 / 2 '0, where a^ = yjfij mw is the harmonic oscillator unit, w is the 
trap frequency, a is the s-wave scattering length of atoms which may be either 
positive or negative, e(t) is the random function on time simulating fluctuations 
of the laser trap frequency, x = il carries the sign of interaction type: "+" 
corresponds to a repulsive interaction (a > 0) and "— " to attractive interaction 
(a < 0), 7 is the coefficient characterizing the damping effects due to influence 
of non-condensed atoms (thermal cloud). 
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Modulation function e(i) can be represented through the deviations of the 
laser field intensity E{t) around its mean value Eq 

<*) = P) 

In order to simplify the model, in what follows we concentrate on the case of 
delta correlated random process with zero mean value: 

(e(t)) = 0, (e(t)e(t')) = 2aH{t-t'), (3) 

where a > is the dispersion of fluctuations, 8(t) is the Dirac delta function, 
and angular brackets stand for statistical average over all realizations of the 
random process e(t). Then the following quantities are adequate to describe the 
evolution of the Bose condensed atomic cloud 

/•OO 

h = 2tt / \4>\ 2 rdr, (4) 
Jo 

Ii = 2tti / (V'Vv — V>Vv) r 2 dr, (5) 



h = 2tt (|V> r | 2 + ||^| 4 )rdr, (6) 

h = 2tt / |^| 2 r 3 dr, (7) 
Jo 

/>OC 

/ 4 = 2ttx / |^| 4 rdr. (8) 



These quantities correspond to the total number of atoms in the condensate, 
radial field momentum, energy of the wave packet, mean square width of the 
atomic cloud and mean field interaction energy respectively. They are governed 
by the system of equations (hereafter a dot stands for the derivative with respect 
to time) 

h = -27/1 + 4/2- 4(1 + e(f))/ 3 , 

h = -27/3 - 2(1 + c(t))Ji - 7J4, (9) 
h = -27/3 + 2/!, 

and the equation for Iq is singled out 

i = -2 7 / . 

The above system is not closed in the presence of the dissipation term, which will 
be considered below. We start with the nondissipative case, 7 = and derive 
the closed system for averaged quantities Jk = (Ik) using for this purpose the 
Furutsu-Novikov formula which gives 

(e(f )/*(*)) = -4a 2 5 M (/ fe (i)) 



3 



(Sk t i is the Kronecker delta). Then the closed linear system of equations for the 
averaged values reads 

j 1= 4J 2 -4J 3 , J 2 = -2J 1 + 8a 2 J 3 , j 3 = 2Ji. (10) 

This system is trivially solved. For the sake of definiteness we concentrate on 
the quantity J 3 which in the case at hand can be interpreted as a mean square 
width of the wave packet. Also, we consider the case of small noise intensity a 2 , 
consistent with the real experimental situations, then the solution has the form 

J 8 (t) = Ce 2it + Cie-t* cos(r]t) + C a e -C * sin(»jt), (11) 

where £ and 77 are introduced as a real and imaginary part of the roots of the 
equation for eigenvalues 

A 3 + 16A — 64er 2 = (12) 

(Aj = — 2£, A2 = £ + if], and A3 = £ — ir/) associated with the linear system 
([To|). In particular for small dispersion a 2 -C 1 we have £ = 2er 2 + 0(cr 4 ), 
7? = 4 + <9(o- 2 ). 

Constants of integration Cj are determined from initial conditions. 

J 3 (0) = x 2 , J 3 (0)=0, J 3 (0) = 8J 2 (0)-8J 3 (0). (13) 
If a Gaussian ansatz is assumed for the initial shape of the BEC 

^(r,t)=A(t)™ P (- 5 ^) l (14) 

where, the amplitude A and width x are connected through normalization con- 
dition for the if) function, we have the following values for constants 

V x V x '/ 

Then it follows from ([ill), that the noise enhances the effects in the BEC. Noise 
induced expansion and contraction of the condensate are described by the same 
eq. (|TT|) , "+" or "-" being assigned to £, correspondingly. These conform with 
positive or negative value of the energy J 2 , therefore the sign in front of the last 
(noise) term in eq.(fl2"|) will be different in these two cases. The noise results 
also in doubling of the frequency of oscillations of the condensate. 

3 Numerical modelling 

In order to verify predictions of the analytical approach we performed direct 
numerical solution of the CP equation (|l|), as well as stochastic ODE's (^), 
resulting from the moments method. Due to increasing role of effective proce- 
dures for integration of nonlinear and stochastic PDE's, we briefly outline the 
numerical approach which has been used. 
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- 2 </>- |0| 2 0-/3 = 0, (15) 
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To proceed with evolution problems at first we find the ground state solution 
to GP equation (Q) with stationary trap potential (e(t) = 0) in the absence of 
the damping term (7 = 0). The corresponding ODE has the form 

d 2 i\) 1 dip 
dr 2 r dr 

where f3 = 2n/tkv is the normalized chemical potential of the condensate. In 
order to solve this equation the following boundary conditions are applied: 
•0(0) = const, ip (0) = 0, which state the regularity and continuity of the ip 
- function at r — 0. Taking into account the fact, that asymptotic behavior of 
the ip - function follows that of the 2D harmonic oscillator solution, one can find 
the value for "0(0) as described in Ref. M. To this end we applied the nonlinear 
equation solving procedure ZEROIN from Ref. || , which combines advantages of 
the bisection and secant techniques. After the both ■0(0) and i[> (0) are known, 
it is strightforward to solve the equation using, for example Runge-Kutta pro- 
cedure. It is worth to mention, that this approach reproduces results for the 
ground state wave functions presented in Refs. @, [lofl . 

The obtained ground state solution has been used as an initial condition 
for modelling of the time dependent GP equation (^) with randomly varying 
trap potential. In order to solve this PDE we employed the method of lines, 
in which the equation is discretized along the spatial variable. From 400 to 
800 grid points are taken to discretize the spatial variable within the interval 
r e [0..5]. Stability of the ground state solution is checked by propagating eq.(0) 
with e(t) = 0, 7 = up to t — 50 in dimensionless time units. No detectable 
variation of the ground state solution is observed. The resulting set of ODE's 
has been solved by the procedure DOPRI8 JTTJ , which is based on the Runge- 
Kutta method of 6(7)th order. This procedure with adaptive step-size control is 
proved to be efficient when providing the given accuracy in a less processor time 
is a crucial requirement, which is the case in solving of stochastic equations. 

The system of equations (|^) when the damping term is neglected, can be 
reduced to the singular Hill equation by introducing the variable X — \/^3 : 

X + 4(l + e(t))X + Q/X 3 = 0, (16) 

where Q = 2I2I3 — I 2 /2 is a constant, determined by integral quantities or the 
initial data. The solution of this equation can be represented in analytical form 
only for particular types of the modulation function e(t). For a stationary trap 
potential (e(t) = 0) the solution is 



X(t) = ^cos 2 (2t) + -J sin 2 (2£). (17) 

In the case of randomly varying trap potential the equatio n (116] ) has to be solved 
numerically. We solved it using the procedure DOPRIN |fll|, which is based on 
the Runge-Kutta method of 7(8)th order. 
In Fig.l the shapes of the trap potential 

U(x) = 2X 2 ~^, (18) 
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corresponding to equation ( |l6| ) at e(t) =0 is shown for cases of repulsive and 
attractive condensates. From this figure it is apparent the behavior of the above 
two types of condensates. In the "particle in a potential well" representation, 
departures of the particle from the stationary point (the minimum of the poten- 
tial well), correspond to expansions and contractions of the atomic cloud with 
positive a, though it remains confined within the trap. Falling of the particle 
into the center of a singular potential corresponds to the rapid convergence of 
the condensate's width to zero (collapse of the condensate with negative a). 

Fig. 2 illustrates expansion of the condensate with a repulsive interaction 
(a > 0) due to random variations of the trap potential. The gain parameter £ 
and oscillation frequency 77 determined from this figure, which is the result of 

Figure 1: The shape of the potential curve eq. ( |Tg| ) for attractive Q = — 1 
(dashed line) and repulsive Q = 1 (solid line) condensates. 

numerical solution of eq.(|l6|), averaged over 400 realizations of random paths, 
are in close agreement with predictions of the moments approach eq. (|Tl|) . The 
important conclusion coming out of this result is the existence of the stochastic 
parametric resonance in the system. Indeed, exponential growth of the am- 
plitude and doubling of the oscillations frequency are clear signatures of the 
stochastic parametric resonance. 



Figure 2: Expansion of the condensate with repulsive interaction caused by 
random variations of the trap frequency. Solid line is according to formula (|ll|) , 
dashed line is the numerical solution of eq. (pf) . X(0) = 0.707, X(0) = 0, a — 
0.1, Q = l. 

The evolution of the condensate wave function under randomly varying trap 
potential, obtained by direct numerical solution of the GP equation ([!]), corre- 
sponding to a particular random path, is shown in Fig. 3. Exponential growth 
of the condensate's width supplemented by doubling of the frequency of oscil- 
lations, is evident also from this figure. 

The time dependence of the mean square width of the condensate, experi- 
encing collapse under the action of randomly varying trap potential, is presented 
in Fig. 4. Calculations are based on numerical solution of eq.([l6|), averaged over 
the ensemble of 400 realizations of random paths. In this case we have the 
interplay between two opposite effects: contraction of the condensate due to 
collapse and its stochastic expansion due to noise. As it is apparent from this 
figure, increasing of the noise intensity a 2 leads to slowing down of the collapse. 
At some level of the noise intensity, contraction of the atomic cloud may be 
balanced, i.e the collapse can be arrested. 
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Figure 3: Evolution of the condensate's profile under the action of random 
variations of trap frequency, a) Numerical solution of eq.(^J) with 7 = 0, \ = +1 
for a particular set of random numbers, b) Variations of the condensate width, 
corresponding to this path, a = 0.5. 

Figure 4: Contraction of a 2D condensate with negative s-wave scattering length 
under different strengths a of the optical trap's noise. Initial conditions: X(0) — 
1,X(0) =0,Q = -1. 

4 Conclusion 

In conclusion, we have studied the dynamics of a Bose-Einstein condensate 
confined in an optical trap with randomly varying trap potential. Applying 
the moments method the stochastic parametric resonance in oscillations of the 
condensate width has been revealed. It is also established, that the noise induced 
expansion of the condensate with attractive interaction between atoms, can 
delay or even arrest the collapse. 
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